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Abstract 

For a given 8 € (— 1, 1), we find out all parameters a, f3 G {0, 1} such 
that, there exists a linear combination of Jacobi polynomials J£+i (x) — 

CJn' (x) which generates a Lobatto (Radau) positive quadrature for- 
mula of degree of exactness 2n + 2 (2n + 1) and contains the point 9 as a 
node. These positive quadratures are very useful in studying problems in 
one-sided polynomial L\ approximation. 
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1 Introduction 

Let a be a positive measure on [—1,1] such that the support of da contains 
an infinite set of points. Fix points —1 < < x n< 2 < ■ ■ ■ < x n<n < 1 and 
consider a quadrature formula of the form 

l n 

f(x)da(x) = Anj/(Xnj) + Rn(f), (1) 
1 3=1 
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with X n _j > 0. For < m < n, this formula is called a positive (2n — 1 — to, n, da) 
quadrature formula [5J, if R n (P) = 0, for all P € P2„_i_ m (as usual, P„ denotes 
the set of all polynomials of degree at most n). The number 2n — 1 — to is the 
degree of exactness. 

We say that a polynomial P n 6 P n generates a positive (2n — 1 — to, n, da) 
quadrature formula, if P n has n simple zeros x n ,i < x n .i < . . . < x n , n in (—1, 1) 
and the interpolatory quadrature formula based on the nodes x n j is a positive 
(2n — 1 — to, n, da) quadrature formula. 

We say that a polynomial Q n +\ generates a Radau positive (2n+l, dx) 
quadrature formula, if Q n +i has n+1 simple zeros y n +i.i < lln+1.2 < • • • < 
2M+i,7i+i in (-1, 1) and 

/l n+1 
P{x)dx = u n+lfi P(-l) + Mn+l,fe-P(2/n+l,fe)) f° r a U P € P2n+1, (2) 
- 1 J 1 



fc=l 



where the weights /Li n +i,fc (0 < fc < n + 1) are positive. For the node 1, we have 
a similar definition. In this case 



/l n+1 
P(x)dx = ^ ^n+l,kP(y n +l,k) + A*n+l,n+l-P(l)) for all P £ P 2jl +1, 
- 1 fe=l 



(3) 



Wc say that a polynomial Q n +i generates a Lobatto positive (2n+2, n+1, dx) 
quadrature formula, if Q n +i has n+1 simple zeros y n +\,i < Vn+1,2 < ■ ■ ■ < 
y n+ i,n+\ in (-1, 1) and 



1 n+1 



P(x)dx = H n+1 . P{-l) + V ^n,kP{Vn,k) + U„, n+2 P(1), (4) 



fc=l 



for all P G P2n+2, where the weights /i n +i,k (0 < fc < n + 2) are positive. 

We remark that, in studying some problems related with best one-sided 
approximation (see [I]), there appear positive quadrature formulas related with 
linear combinations of Jacobi polynomials of the form 

T^\x) = J^\x)-CJW){x), (5) 

where jjf'^\x) is the monic Jacobi polynomial of degree k. In particular, when 
we want to characterize the polynomials of the best one-sided Li-approximation 
to Heavisidc functions (in this case 9 6 (—1,1) is the parameter defining the 
Heaviside function), the following problem arises: Does the polynomial P„+f 
defined by (J5j generate a Lobatto (Radau) type positive quadrature formula? As 
we will show here the analysis is not simple, because for n > 1 and a, G {0,1} 
fixed, we can have such a positive quadrature only in some specific subintervals 
of (—1, 1), while we need to cover all the open interval (—1,1). Therefore, some 
changes of the parameters a and (3 are necessary, and we need to know when and 
how the changes should be done. These facts were the motivation for writing 
this paper. 
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Positive quadrature formulas associated to linear combination of orthogonal 
polynomials have been studied by different authors (for instance, see [6] and 
[7J and the references therein). In particular Theorem 4.1 of [7J provides suffi- 
cient (and complicated) conditions to obtain Lobatto and Radau type positive 
quadrature formula associated to Jacobi weights, with a,/} £ {0, 1}. In [7] the 
author deals with general orthogonal polynomials. Here we need a more careful 
analysis. That is the reason why we pay attention to Theorem 2.1 of [7J, instead 
of Theorem 4.1. The use of special properties of Jacobi polynomials helps us to 
obtain more specific results. Moreover, we will present the result it terms which 
depends on the parameters a and (3 defining the Jacobi polynomials. 

The paper is organized as follows. In Section 2 we collect several known facts 
related with Jacobi polynomials which will be needed. In Section 3, for n G No 
and a,/3 S {0, 1}, we first find all of real C for which the polynomial T^"'f ^ given 
in (J5J generates a positive (2n,n + 1, (1 — x) a (l + x) l3 dx) quadrature formula. 
Some of these C generates a Lobatto (Radau) type positive quadrature formula, 
others not. For the case when a + j3 ^ 0, wc find out all C for which there is 
a Lobatto (Radau) formula. In the last section we consider a more complicated 
problem. For 8 G (—1, 1) and n G No fixed, we want to find out all Lobatto 
(Radau) positive quadrature, with degree of exactness 2n + 2 (2n +1), which 
contains 9 as a node. Of course, the problem has interest only when 8 itself is 
not a zero of a Jacobi polynomial of degree n (with a,f3 G {0, 1}). That is the 
reason why we define 

^ ) W = 4Sl7^ M^O^O), (6) 

Jn (X) 

and consider the linear combination 

= J&?\x) - Qi a ' 0) (e)Ji^(x), (7) 

where 8 G (—1,1). The question is: for which selections of a and f3 (a, (3 G 
{0,1}) the polynomial R^.^g{x) generates a Lobatto (Radau) type positive 
quadrature? This is just the problem which appears in studying best one-sided 
polynomial Li-approximation for the Heavisidc functions. Wc will solve it. For 
the solution we use some recent results related with the interlacing properties 
of zeros of the Jacobi orthogonal polynomials given in [5] . 

2 Some facts related to Jacobi polynomials 

For a > — 1,/3>— 1 and n G N, the Jacobi polynomial Pi a '^ is the unique 
polynomial of degree n which satisfies 

J Q n ^{x)P 7 [ a ^(x)(l-x) a {l + x^dx = 0, forallQ„_i eP„_i, (8) 

and 

pi-«d) = ( n ; a ). o) 
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It is known that (see [5], p. 59) 

P^\x) = (-i r p(M(- x ) and P^\-l) = (-1)" (^ + J 

The coefficient of x n in Pn"''^ (see [5], p. 63) is given by 



(10) 



2" 

The monic Jacobi polynomial will be denoted by jjC'^\x). Since 



4 a ' P \*) = -r^Pli a ' P) W> (12) 

in 

one has 



T (~,P)m= 2 T T(n + a + /3 + l) r(n + a + l) 
" lj r(27i + a + /? + l) r(a + l) [ ' 

and 

r(q,ff)/ n _/ n „ 2^T(n + a+g+l) g(n+g+l) ^ 
J " 1 iJ r(2n + a + /3 + l) r(/? + i) 1 J 

The Jacobi polynomials satisfy a recurrence relation (see [9], p. 71): 
P (Q ' /3) (x) = 1, P 1 (Q ^ ) (x) = [{a + + 2)x + a - 0\/2 
and, for n > 2, 

2n(n + a + /3)(2n + a + /3-2)P^(x) (15) 

= (2n + a + /3 - l)[(2n + a + /3)(2n + a + /3 - 2)x + a 2 - /3 2 ]P^"f ' (x) 

-2(n + o- l)(n + /3- l)(2n + a + ^)P^(x). 

For the monic Jacobi polynomials Jn*'^\x), with n > 2, one has (see [2], p. 
153) 

4-«(x) = (x - c^>) j£f >(x) - X^J^\x), (16) 

where 

c («>/3) = ^ ~ a ' (17) 

(2n + a + /3-2)(2n + Q + /3) V 1 

and 

X ( a ,0) = 4(n - l)(n + a - l)(n + - l)(n + a + /j - 1) 

(27i + a + /3-2) 2 (2ra + a< + /3-l)(2ra + a + /3-3)' 1 j 

When n = 1 we should take c^' p) = (j3 - a)/(a + j3 + 2). 

From HI]), it follows that = ln' a) ■ Hence, from fflOl) and flUJl one has 



!n 1 ' ~ p(a,/J), * ~ ,(a,/3) p08,a), x ~ ^» W' ^ 
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Moreover, from (fT5|) and (HU), we obtain 



><<*.« m = 2(n + a + 1)(n + a + /3 + 1) > o f2o) 

W {2n + a + /3 + 2)(2n + a + f3 + l) { ' 



and 



q^(-d=-^(i) = - j:yj?tti^+B l l u < °- <™ 

(2n + a + p + 2){2n + a + p + 1J 
A direct calculation shows that 

It is known that the zeros of the Jacobi polynomials are real, distinct and 
are located on the interior of the interval [—1,1] (see [5], p. 44). For n £ N, the 
zeros of the Jacobi polynomial Pn^ will be denoted by x„ jf^ (1 < k < n), 
with 

-Ka^<-<»fe«<l- 
The zeros have the interlacing property (see (S], p. 46 or [2], p. 28): 

_(<*.« ^ _(<*.« <- r («>« ^ -(«,« <• T (a,« ^ _(«,« r2 q\ 

x n+l,l ^ ■''n,! v x n+l,2 - x n,2 ^ ■ ■ • ^ ^ •*'n+l,tt+l" y^ J 

In what follows, for simplicity we use general parameters a., ft > —1, but we 
need only the cases a, (3 £ {0, 1}. 

3 Quadrature formulas 

Proposition 1. For a,/3> —1, n £ N, the polynomial T^^\x) given in ([5]) 
generates a positive (2n, n + 1, (1 — x) Q (l + x)@dx) quadrature formula if and 
only if 

Q^X-l) < C < Qt?Xl). (24) 

Proof. It follows from c) in Theorem 2.1 of [7] (with m = 1) that T^'f^x) 
generates a positive (2n,n + 1, (1 — x) Q (l + x) l3 dx) quadrature formula if and 
only if the following conditions hold 

i) there exists a constant K, with \K\ < 1, such that 

T^\x) = (x + K)J^\x) - Jtf >(*), 

where A^+f' is given by (fT%|). and 

ii) sgn(T I ^f ) (±l)) = (±l)«+ 1 . 
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Notice that, if J^ a,0) (x) ^ 0, then 

= J&f (*) -CJ^\x) = J^\x){Q^\x) - c). 

So, 
and 

T^f ) (-l) = J^)(-l)(Q^)(-l)-C). 

From (|14[) we conclude that the condition ii) is just equivalent to (|24[) . So, 
we only need to prove that ([24)) implies i). 
Indeed, from (|I6j) . we can write 



-^n+'f ^ (* c ) — ( x c n+l ^)^n { x ) -^i+l ' J n —\ ^ i x ) ■ 

Thus, it is sufficient to prove that | c^!f + C |< 1. But, from ||22) we obtain 

,(«,/?) 

V n + 1 

(-1) 



-i < -i - g<-«(-i) - ^ - + o<r-«(-i) = + o^»(-i) 



< «&« + c < + gfi».«(i) = i - g<-«(i) - -0- + Q^\i) < i. 

□ 

For our next result we need the following Radau quadrature (see 0J, p. 169). 
Theorem 1. For each in G N and every polynomial P € P2m, o^e /ias 



/i m 
P(x)dx = A m P(-l) + ^i? m , fe P(^), 
1 k=l 

/I m 
P(x)cfe = P m P(l) + ^S m , fe P 
" 1 I 1 



fc=l 

and 

/•l m 

.(1.0) 

fe=l 

where, for 1 < k < m, R mi k and S m ^k are positives and 

2 



X m,k )i 



(25) 



(26) 



{m + 1) 2 ' 

Theorem 2. Pis; neE 

i) The polynomial t'^V (x) generates a positive (2n, n+1, (l + x)dx) quadra- 
ture formula and a Radau positive (2n + 1, n + 1, dx) quadrature formula 
as in ([2]) z/ and onZ?/ if 

*"<- 1 ) < -s^ < c < " IS' (27 > 
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ii) The polynomial (x) generates a positive (2n, n+ 1, (1 — x)dx) quadra- 
ture formula and a Radau positive (2n + 1, n + 1, da;) quadrature formula 
as in (j3]) if and only if 

* 0> (- 1 ) = -|S <o< l^ < «» , ° , < 1 >- < 28) 

Proof. We will prove i). The proof of ii) is similar, 
i) (=>) From Proposition Q] and pip wc know that 

- (2„ ( : 3 + )(f +1) = « "f- 1 ' < c < - is- « 

Moreover, it follows from ^ that 

Taking into account that (see Theorem [I}, for each m G N, 

(m + 1) 
one has 

I ^ (^ « M) - C VTW 40,1)( - 1} 

2jf 1) (-l) / (n + l) 2 (oa) r \ 2jf'(-l) / n+1 ' 

(n + 1) 2 ^(,i + 2) 2 W " 1 } °y (n + 1) 2 V 2n + 3 
On the other hand 

= jft?(-l)-CjS , ' 1 )(-l)= J^-lXflg^-lJ-C). 

Therefore 

1 ( n+1 _ 

(n + l)#(-l)-C)l 2n + 3 

Since Q { n ,1} {-1) - C < and ^ n+l! o > 0, we conclude that 

n+1 



C > 



2n + 3' 



(<=) If ([27)) holds, then the polynomial T^V^ generates a positive (2n, n + 
1, (1 + x)dx) quadrature formula for n > 1 (Proposition [T]) . Hence, for every 
polynomial P £ P2 n 

/l n+1 
p(x-)(i + x)dx = K+i,kP{y n +i,k), (30) 
fe=i 
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where the points {j/n+i,j}^Ji are the zeros of T$_{> and 

Xn+x,k = / h(x)(l + x)dx, 



where £k(x) are the cardinal Lagrange polynomials at the nodes {j/n+i,j}™=i- 
Moreover, from (|30]) we have that, for each S n -i £ P rl -i, 

1 

S„_ 1 (z)T^ ) (.t)(1+.t)o':e = 0. (31) 

-1 

Note that, every polynomial P £ P n +i can be written in the form 

^,(0,1)/ s n+l p, n 

Hence, from ([30|) . for each P £ P n +i, we have 

/l n+l 
P(x)dx = j Lt n+ i )0 P(-l) + Mn+l,j'P(yTi+l,*)j 
- 1 fe=l 

where, for 1 < j ' < n + 1, 

fJ-n+l.j = 1 . > and Mn+1,0 = r ( ,i^ ~ J i T r ( l +l ) ( a; ) da; > °> 

where for the last inequality we have used the computations given in the first 
part of the proof. 

Finally, for every Q2n+i £ P2n+ii there exist polynomials 5,1-1 £ P«— l and 
Pn+l £ Pn+1 such that Q 2n +i(x) = S n - 1 (x)(l + x)T^_\ ) (x)+P n+1 (x). Taking 
into account pip , wc obtain 

i ,i 
Q 2n +i(x)dx = I P n+ i(x)dx 
-l J -i 

n+l 

= (J*n+lfiPn+l(— 1) + Mn+l.jPn+l (g/n+l.fc) 
fc=l 
n+l 

= Mn+l,oQ2n+l(— 1) + A i n+l,jQ2n+l(yn+l,A;)- 
fe=l 

□ 

For the next result see [1], p. 172. 
Theorem 3. For m £ N and every polynomial P £ P2 m +i, one ftas 

/l m 
P(x)dx = ^ A m>fc P(^;p) + S ro P(-l) + C m P(l), (32) 
_1 fc=i 
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where, for 1 < k < m, A m ^ are positives and 

2 



(m + l)(m + 2) ' 



Theorem 4. For n 6 N, £/ie polynomial T^Vp(x) generates a positive (2n,n + 
1,(1 — x 2 )dx) quadrature formula and a Lobatto positive (2n + 2,n + 1, da;) 
quadrature formula as in (01 i/ and only if 

" + 1 <C<^. (33) 



2n + 3 2n + 3 

Proof. ( =>) From Proposition [1] we know that 



g(M) ( _i )<C<g <M) (1) = £±3 (34) 



2n + 3 ^ n y ' ^ n y ' 2n + 3 

Moreover, it follows from ((4]) that 



1 



Mn+i,o = (i 7T- / T^ 1 ) (x)(l -x)dx > 



and 



1 



^n+l,n+2 — 



T^(x)(l + x)dx > 0. 



L n+i 

From Theorem [31 for each m £N, we have 



and 



Then 



4^(^(1 + ^= (m + 1)(m + 2) ^(D, 



1 j£< 1 \x)(i-x)dx= 4 ja.D ( -i). 

.! (m + l)(m + 2) 



T&ftsXl + x)dx = 1 1 (^(x) - CJ^x)) (1 + x)dx 

= 4 j^m - 4C jWm 

(n + 2)(n + 3) ™ +1 1 ; (n+l)(n + 2) " 1 J 

4ji M) (l) / n+1 o(1 , 1)n s _ r \ 4J' M )(1) / n + 1 _ > 
(n+l)(n + 2) \n + 3 V " lj 7 (n + l)(n + 2) V2n + 3 , 

and, similarly, 



(n+ l)(n + 2) V 2n + 3 
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Hence 

4 ( n+1 _ 

(n+l)(» + 2)(g^(-l)-0) I 2n + 3 

and 

4 f n + 1 _ 

M " +1 ' rl+2 ~(n + l)(n + 2)(Q^(l)-C) U + 3 ° 

Taking into account (|34[) we obtain that /i„+i,o > and /i„+i, n +2 > if and 
only if (021) holds. 

(<=) If (pJ3"|) holds, then, in particular, (jM)) holds and then, from Proposition 
[U the polynomial T^T^ generates a positive (2n, n + 1, (1 — x 2 )dx) quadrature 
formula for n > 1. Hence, for every polynomial P € P2„ 



/l n+l 
P(x)(l - x 2 )dx = ^ A n+ i jfc P(y n+ i jfc; ), 
1 fc=i 

where the points {yn+i,j}?=i are the zeros of T^fi and A„, fc > are given by 
A„+i,fc = / £k{x){l - x 2 ) dx, ft = 1, ...,n+l, 



where are the cardinal Lagrange polynomials at the nodes 

Moreover, for each S € P n -i, 



l 

S{x)T^{x){\ - x 2 )dx = 0. (35) 

-l 

Note that, every polynomial P 6 P„+2 can be written as 

2T n+l (- 1 ) fc=l 

+P(1) d + ^:'w . 

2r^;>(i) 

Hence, for each P e P n +2 , 

/l n+l 
P(x)dx = (j, n+lfi P(-l) + n n+ x >k P(y n+ i >k ) + / u„ + i jn+2 P(l), 
- 1 fe=i 



where, for 1 < ft < n + 1, 



A n +1 



Mn+l.fe — ^ T > U ' 

1 - Vn+l.k 
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1 f 1 

Mn+i.n+2 = — / ri+r(a:)(l + a:)da; > 

(1) 

and 

1 f 1 

where for the last inequalities we have used the computations given in the first 
part of the proof. 

Finally, if Q2n+2 G ^2n+2, then there exist polynomials S n -\ G P n _ i and 
Pn+2 G P n +2 such that Q 2 „ +2 (:r) = S n _i(x)(l - x 2 )T^\\x) + P n+2 (x). Taking 
into account (1351) we obtain 



Q2n+2(a;) c?a; = / Pn+2(x)dx 

l J-i 

n+l 



= Mn+l,0-Pn+2( — 1) + J]] Mn+l,fc-fri+2(?/n+l,fc) + Mn+l,n+2-Pn+2 (1) 
fc=l 
n+l 

= A t n+l,oQ2n+2( — 1) + /J Mn+l,fc £?2n+2 (Vn+IM ) + A*n+l,n+2<92n+l (1) • 



fc=l 

□ 

4 Positive Lobatto and Radau quadratures con- 
taining a give node 

For a,/3 > -1, let {xffi}^ be the zeros of P& a ' P \x). Put cc^o /3) = -oo and 
x ^nn+i = +°°- It is know that the function Q^*'^\x) defined in © increases 
from — oo to +00 in each one of the intervals (x^f_} 1 ,x^f > ), 1 < k < n + 1, 
(see Theorem 3.3.5 of p. 46), and Qn is negative on (x^nk-ii x ^n+ik) an ^ 
positive on ■ 

Taking into account the above properties, we define some numbers as follows. 

For each k, 1 < k < n + 1, let i>„+i fc e ( x nk-i! fc) be the om y P om t 
where 

Q^ ) {^ >k ) = Q^ ) {-i) (36) 

and let u^+f * G (^i+i^j^nfc ) be the only point where 

a ( ^ ) K ( l+ fi) = Qi Q ' /}) (i)- (37) 

, (Q ' /3) - -1 and 7/, (Q ' /3) 

J ri+M — 1 ana W n+l,n+l 

We define the intervals -Z^+'f j. = ( u n +i fe) fe ) , = 1, . . . , n + 1. Since 

(a,/3) (a,/J) (a.fl , _ fl 

X n,fc-1 < V n+l,k < X n+l,fc ^ W n+1M ^ X n,fc ' ™ U, . . . , 7i, 
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the function R£$ (x) (given by ©) is well denned for 9 6 {Jktl 4.+i,fc- Morc - 
over, since Q n is increasing on (x„ k _ 1 ,x n k '), we obtain 

Hence, we can rewrite Proposition [T] in terms of the sets /^f^. 

Proposition 2. Fix a,j3> —1, n E N, 9 E (—1, 1) and 1 < k < n. The polyno- 
mial Rn+i\( x ) defined by generates a positive (2n,n+l, (1 — x) a {\ + x)@dx) 
quadrature formula which contains 9 as a node if and only if 6 € Ufe=o ^i+ffe- 

Notice that, for a fixed a, j3 6 {0, 1}, we can not cover all the interval (—1, 1) 
with the union of the intervals considered in Proposition [2] Hence, if we want 
to find a positive quadrature formula for every 9, the parameters a and /? must 
change. 

We will prove that the union of the sets I^l with a, /3 € {0,1} covers 
the interval (—1,1). When 9 belongs to one of these intervals, we will find a 
positive quadrature formula associated to a linear combination of some Jacobi 
polynomials, which contains 9 as a node (see Proposition From Theorems 
[2] and [3J we know that the interval where we can have a Lobatto or a Radau 
positive quadrature formula is a proper subset of an interval where we have a 
quadrature associated with a Jacobi weight. Thus we should verify that some 
of these intervals have non empty intersection. 

Taking into account Theorems [D and [3J for a, (3 E {0, 1} and 1 < k < n + 1, 
we also consider the points defined by the equations 

2n + 3 ' 



Vn+l,k) n _i_ o 1 Vn+l,k ^ ( U n+l,fc' X n+l,k I ' 



V"n+l,k) ~ 2n + 3 ' n + 1 ' k \ w n+l,kJ ■ 

From a Theorem of A. Markov (see Theorem 3.3.4 of [5], p. Ill) we know 
that the function 

F„, fc (a,/3) =4?^' l<fc<n, a > -1, /3 > -1, 

decreases with a and increases with (3. Moreover, from Theorem 2.3 in [3] we 
have, for 1 < k < n, 



(a,/3) ^ ^(a+s,/3+t) ^ ^(a 
n+l,fc 



< < k < <+i fc+i> for all a, > -1 and < s, i < 2. (38) 



We do not know if the interlacing property given in (|39[) is a new result. 
Proposition 3. If n E N, then 

i) 



Vnfh+l = X n'k = W n+lk > 1 < fc < »> 
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and 

W n+lk = X n;k and v £3,W= X ni ) > I < k < n. 

In particular, 

X n,l < X n,l < X n,2 < X n,2 < ■ ■ ■ < x n,'n) < X n,'n) • (^9) 

nj For 1 < fc < n + 1, one has 

T (i.o) . (0,0) _ (1,0) (1,0) (1,0) 

•''n+l.fc ^ x n+l,fc ~~ ^n+l,k v "Wl,fc ^ x n,fc ' 

(0,1) (0,1) (0,1) _ (0,0) (0,1) 

x n,fc-l ^ u n+l,k ^ IJn+l,k ~ x n+l,k x n+l,k ' 

r (l,l) . (1,1) <- _ T (i.o) . (1,1) 

x n,k-l v ^n+ljfc v Wn+l,fc — x n+l,k ^ x n+l,fc ' 



and 



T (i.i) ^ r (o,i) _ Ji.i) . ,,,(i,i) . 

■^n+l.ft ^ x n+l,fc — ^ u 'n+l,fe x n,k 



Proof, i) It is known that (see equation (4.5.4) in [5] p. 71) 

,_ 2(n + l) P^\x) - Pj° + \\x) 
" W ~ 2n + 3 1-z 



t(0,1)/_(1,1)x ,(l,0) p (0,l) r (1,1)x 
J ™ V x n,fc ) l n+l r n V-'n^h I 

□ (O.lVlN 
n+1 W _ /-)(0,1)/ 



,(o,i) ,(o,i) p (o,i) rn 
,(o,i) 7 (o,i) D (o,i)/— Wn 1 j - 



(40) 



and 

p a,» (x y_ *(n + l)P£> \x)+P™(x) 
" 1 j " 2n + 3 1 + z ' 1 J 

In particular, from (|4"0")) and (|4"T1) we have 

-fn+l^( x i,fc^) — -Pn ' 1 ^^!,^;^) an d ^n+1 (^n.fc ~ — ^n 1 '°H a 'i,fc 

Therefore, for 1 < fc < n, one has (sec ([9])) 

r (o,i) r (l,l)s 7 (0,l)p(0,l)/ (1,1)-, 
^(0,1)/ (!>1)\ _ J n+1 \ x n,k ) _ Ln r n+l \ x n,k J 
' \ X n,k ) 



L n+1 l n+l r n \ L ) 

Since the zeros of Pn~'^ are symmetric, we also have Qn (—x^'P) = Qn :1 \l). 
Then, from (fT9]). we also obtain 

Moreover, it follows from the first inequality in (f3"5)) (with a = 0, P = 1, s = 1 
and < = 0) and Markov's theorem that 

»(o,i) ^ „(i,i) ^ _(o,i) 
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Since Qn' 1 ^ is increasing in the interval k-> x n'k^)> we 

W n+l,k = X n,'k ' k = l,...,n. 

On the other hand, from Markov's theorem and the second inequality in ()38l) 
with a = 1, f3 = 0, s = and t = 1, we know that 

T (i,o) ^ (i,D ^ T (i,o) 
Since Qn is increasing, we have 

V n+l,k+l = X i,fc i & = 1> ' ■ • ) n - 

It follows from equation (4.5.3) in [9] (p. 71, with a = /3 = 0) that 

p(0)0) / \ p(0,0) , , p (0,0) , x p (0,0) / s 

1 — a; 1 + x 

In particular, from (|42[) we have 

Therefore, 

.(0,0)/ (1,0) x ,(0,0) p (0,0), (1,0) v (0,0) 

0(0,0)^(^0) n _ J n+1 ( X n,fc j _ 6 » fVj-1 ^n.fc J _ In -n( -°)(l\ 

Hn K ^ .(0,0) , (1,0). ,(0,0)p(0,0)/ (l,0)s ,(0,0) ^« 

•'n ^ x n,fc / 'ra+l \ x n,h > l n+l 

Moreover, from the first inequality in (p?5|) (with a = f3 = 0, s = 1 and t = 0) 
and Markov's theorem, one has 

T (o,o) , (1,0) (0,0) 

and then w£$, k = a££\ k = l,...,n. 
In a similar way, we obtain 

From Markov's theorem and the second inequality in (|38|) (with a = 0, /3 = 0, 
s = and t = 1) we know that 

T (o,o) (0,1) (0,0) 
^n,k ^ x n,fc ^ x n+l,/c+l' 

and then = 4°fe ) > & = 1, • • • .n. 

Finally, since 

T (0^0) (0,0) _ (1,0) (0,0) (0,0) _ (0,1) (0,0) 

x n+l,fc ^ W n+l,k ~~ X n,fe ^ x n,fc ^ "n+l.fc+1 — X ri,fc ^ x n+l,fc+l' 
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we have proved 

ii) From the recurrence relation (1151) we obtain 



p(0,0), (0,0) \ 
r n+2 \ X n+l,k) _ n+1 



Thus (HU) and (gSJ yields 



.(1,0), (0,0) x (i i0 ) p (i,o) r (0,0) x 

^(1,0)/ (0,0) \ J n+1 y^n+l.k) _ <m r n+l \ x n +l,k> 

V„ { X n+l,k) r(L0)/ (0,0) x ~ ,(1,0) p (l,0)/ (0,0) x 

•Jn \ x n +l,k) l n+l r ™ \ x n+l,k) 



, p(0,0)/ (0,0) x _ p (0,0), (0,0) v 
71 + 2 Jjt+l ( X ,t+1,/J \ x n+l,k) 
2n + 3 p(0,0), (0,0) v _ p (0,0) f (0,0) v 

, o p(o,o)/ (o.o) , 

^ + 2 "^+2 l^n+l.fc) "+1 



2n + 3 p(o,o)/ (o,o) , 2n + 3 



= Qi 1 '° ) (4?)- (43) 



From Markov's theorem and the first inequality in (|3"51) (with a = (3 = 0, 
s = 1 and t = 0) we know that 



T (i.o) < -(0,0) . (1,0) 

x n+l,fc ^ x n+l,fc •''n,*: ■ 



Since Qn increases in {x^+lk-> x^'^), one has 



(1,0) _ (0,0) (1,0) ,, _ 1 j 

Since the zeros of J 5 ^? are symmetric and the numbers in (|43[) docs not 
depend on k, we also have 

n(i.o)f_ T (o,o) ] _ n± 1 
I ^+1,J - 2n + 3' 

Then, from (fl9|) we obtain 



( n + 2 ) 2 = n(o,i) ( '_i') = 0(0,1)^(0,1) v 

(n + l)(2n + 3) 1 J W " 1 " +1 < fcJ 



From the second inequality in (|38|) (with a = /3 = 0, s = and t = 1) and 
Markov's theorem we have 

-(0,1) < (0,0) (0,1) 
x n,fc-l ^ x n+l,fc ^ *«+!,/;• 



Hence, yl l+ i k — x ^ n+ i^ € (^nli^j^n+i.fc)' 
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From (|4ip and the recurrence relation (fT5j) we know that 

p(V)/„(l,0) s , , _ u - , , p (i,o) r (1,0) n , p (i,o) r (1,0) x 
fn+l \ X n+l,k> _ [n + 2)(2n + 3) Vi gn+^fcj + \ X n+l,k> 

p(i.i)/ (1,0) x ~ (n + l)(2n + 5) p( 1 '°)( T ( 1 .°) \ + P (i.o)/ (i,o) x 

(n + 2)(2n + 3)^ ) (4 1 ; i, fc )) 
(n + l)(2n + 5) P a.o) (s (i,o) jb)) 

(n + 2)(2n + 3) (n + l)(2n + 5) _ n + 2 
~~ ~ (n + l)(2n + 5) (n + 3)(2ra + 3) ~ ~n + 3' 



Therefore 



,(1,1) p(l,l)/ (1,0) X 
tra ^n+l \- L n+l,k) 

7(Ml p(M), (1 

(n+l)(n + 3) n + 2 
(n + 2)(2n + 3) ?i + 3 



^„+i^„ + i, fc ; - (M) (i,i) (i )0 ) . 

l ri+l r n l x n+l,fc/ 



1 



~~ ?n + 3~ ^n+i,fe' 1 2n + 3 ~ _ s u n+i,)tJ 

From the second inequality in (|38[) (with a = 1, /3 = 0, s = and i = 1) and 
Markov's theorem we know that 

T (M) (1,0) (1,1) 

Hence, yj^. = sjV? fc , fe = 1, . . . , n + 1. 



Finally, from (|40| and the recurrence relation ([15]) we obtain 

p (i,i) r (0,1) x . p (o,l) r (0,1) x p (0,i) r (0,1) x 

x^n+i.fci (n + 2)(2n + 3) J^+i l^„ + i, fc j ~ ^ n +2 \ x n +i,k) 

p(W)f T M l fn + l)(2n + 5) p(°.i)( T (°.i) \ P (o,i)/_(o,i) x. 

r « x^n+l.fci V ;V ' r n \ X n+l,k) r n+l\ X n+l,k) 

(n + 2)(2n + 3)^ ) (4+ 1 i ) fc )) 
(n + l)(2n + 5) P ^) (l (o,i) fc)) 

(n + 2)(2n + 3) (n + l)(2n + 5) _ n + 2 
(n + l)(2n + 5) (n + 3)(2ra + 3) _ n + 3 ' 

,X) Pn+i\ x n+lk) _ (n + l)(n + 3) n + 2 



Therefore (see (|TTt ) 

(i,i) r (0,1) X 



2n + 3 v 2n + 3 

From Markov's theorem and the first inequality in (|38l) (with a = 0, /3 = 1 , 
s = 1 and t = 0) we know that 

r (i,i) < r (o,i) < T (i,i) 

x n+l,k ^ ^n+l.fc v '''n^ 1 

Hence, = 4+i,fc ,fc=l,...,n + l. □ 
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5 Main results 

Theorem 5. Fix 6 G (—1,1) and n G N. The polynomial R^+i e( x ) generates 
a Lobatto positive (2n + 2, n+ 1, dx) quadrature formula as in Q if and only if 

U \ n+l,fc> X n+lJ 
k=l 

Proof. From Theorem [H we know that the polynomial R^+i g(x) generates a 
Lobatto positive (2n + 2, n + 1, dx) quadrature formula if and only if 

= < o^w < £± = e 1 ' (*£&)■ 

Since Qn increases, it happen if and only if there exists k (1 < fc < n + 1) 

such that P f« (M) z (1,1) \ Rut 7/ (1,1) - r (1 ' 0) and z (1,1) - t (0>1) n 
sucn tnat f g (,2/ n+1)fc) z n+1>k ). am, y n+l k — x n+l k ana z„ +l fc - x n +i,k- u 

Theorem 6. Fix G (— 1, 1) and n G N. TTie polynomial R^+i e( x ) generates 
a Radau positive (2n + 1, n + 1, da;) quadrature formula as in ([2]) «/ and onZj/ i/ 



^(-LxfiMU U 



,(0.0) ^ I Ifl I (1,1) (0,0) 

\^n,fc ' "''ra+l.fc+l 

Vfc=l 

Proof. From Theorem [5J we know that the polynomial R^f^ (x) generates a 
Lobatto positive (2n + 1, n + 1 , dx) quadrature formula if and only if 

Q^i^L) = e o) (-D = -|^| < Ql 1 ' ^) < £^ = Q™(*£&)- 

Since Qn'°^ increases, it happen if and only if there exists k (0 < k < n + 1) 
such that 8 G (t^+i fcs 2 i+i k)- That is, 

n+l / n \ 

6e I I fw (1 ' 0) z (1 ' 0) )-fv (1 ' 0) z (1 ' 0) )l I I I (v ih0) z {1 ' 0) ) I 

" fc 1J V V n+l,fc' Z n+l,fcJ — l W n+l,l' Z ri+l,lJ U I U ^ ra +l,fc+l ' z „+l,fc+l/ I ■ 
k=l \k=l ) 



(!.0) _ 1 ..(b 1 ) _ Jbl) „„j ,(1,0) _ (0,0) 



Rut 7i v ' ' — -1 v y ' ' - x y ' ' anrl 

But, v n+11 - 1, v n+uk+1 - x nM ana 



"n+l,fc — ^n+l.fc' 



□ 



Theorem 7. Fix G (—1,1) and n G N. TTie polynomial R^+\ g ( x ) generates 
a Radau positive (2n + 1, n + 1, dx) quadrature formula as in ([3]) i/ and on/?/ if 



'e (U(*&V£iP) JUG 



(0,0) 
n+l. n+l' 



Vfc=l 



17 



Proof. From Theorem [2] we know that the polynomial R^+i g(x) generates a 
Lobatto positive (2n + 1, n + 1, dx) quadrature formula if and only if 

o^foK*) = < e 1 ^) < = o^d) = Q^SS, k ). 

Since Qn increases, it happen if and only if there exists k (1 < k < n + 1) 
such that 6 £ (y n °+l k , w n+lk)- Tnat is ' 

n+l / n \ 

61 G I I f« (0,1) «7 (0,1) 1 - I I I (v {0 - 1} w {0 ' 1] Mil Cy^'V w (0ll) 1 

fc=l \fc=l / 

But, w { °+l n+l = 1, to^ 1 ^ = x^f, k = l,...n, and y^ 1 ^ = sJJ ^, fc = 
l,...,n + l. 

□ 

Theorem 8. Fix 6 £ (— 1, 1) and n £ N. XTie polynomial R^+i $ ( x ) generates 
a positive (2n, n + l, dx) quadrature formula as in (JlJ if and only if 



* e (-i.* 0) ) U (U U (^ 1J ,i) 



Proof. From Proposition [TJ we know that the polynomial R^+i g( x ) generates 
a positive (2n, n + l, dx) quadrature formula if and only if 

g£' 0) («£S*) = Qi°' 0) (-i) < Ql°' 0) W < Ql°' 0) (i) = Ql^H^V- 

Since Qn' '* increases, the conditions given above holds if and only if there exists 
k (1 < k < n + 1) such that 9 £ (u^ ^, w£ + °i fc )- That is, 



o s (e^.-iri) u (u u (« 



(0,0) (0,0) 
n+l, n+l' w n+l,n+l 



b^, a = -i, = ^ a < * < »). - * o) (!<*<«) 

and = 1. □ 

Remark 1 . A direct computation shows that Theorems [5j [6] and [7] also holds 
for n = 0. 

Corollary 1. For all m £ N and /or a// 6> £ (—1,1) i/iere exisi A: £ No and 

a,/3 £ {0,1} (uniquely determined) such that R^f^ 1 generates a quadrature 
formula of order of exactness m. 
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Proof. Assume m = 2n + 1 with n e N . If = then it is know that 

P^+i generates a Gauss quadrature formula of order of exactness 2n + 1. Since 
Pl°' 0) (6) ^0, we have 

-Ri+l i 6/( ;E ) = Cn,e-P n +i (^), 

for some C ni e 7^ 0, and the result holds by taking a = f3 = and fc = n + 1. 
A similar result follows if is a zero of ■ 

Note that the intervals in Theorems [6] and [7] are disjoints and cover the 
interval (—1,1) unless the points {Xnk^}k=i an d { x n+i k}k=i- So, setting 
x„o = — 1 and a^^Li = 1, we have the following situations: 

(i) e (»n,lb } > 4+°i,fc+i) for somc * = 0, 1, • • • 

(ii) 6> e (^i+i^fcj^i^J f° r somc = 1, • • ■ ,n + 1, 

(iii) = x„ z. for somc k = 1, . . . , n, or 

(iv) 6* = x^V^ fe for some fe = 1, . . . , n + 1. 

In the case (i) Theorem[6]assert us that the polynomial R„+i g provides a Radau 
quadrature formula of order of exactness 2n + 1. 

The case (ii) follows analogously from Theorem [7] by using the polynomial 
Rn+ie- W e also have a Radau quadrature formula of order of exactness 2n + 1. 
In the case (iii), Theorem [3] provides a Lobato quadrature of order of exactness 
2n + 1 that contains the point 9 as a node. This formula has n interior nodes 
and two extremal nodes. 

Finally, in the case (iv) it is know that the zeros of the polynomial P^+f^ 
provide a Gauss quadrature formula of order of exactness 2n + 1. 

The proof follows analogously when m — 2n, n > 1, by using Theorem [5] 
(with n replaced by n — 1) and Theorem [8l when 9 is not a zero of P^ ' 1 ^ nor 
Pn , and from Theorem [1] otherwise. □ 
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